Coulomb dissociation of and the low-energy cross section of the '^Be(p,7)^B solar 

fusion reaction 
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An exclusive measurement of the Coulomb breakup of *B into '^Be+p at 254 A MeV allowed 
to study the angular correlations of the breakup particles. These correlations demonstrate clearly 
that El multipolarity dominates and that E2 multipolarity can be neglected. By using a simple 
single-particle model for *B and treating the breakup in first-order perturbation theory, we extract 
a zero-energy S factor of 5i7(0) — 18.6 ± 1.2 ± 1.0 eV b, where the first error is experimental and 
the second one reflects the theoretical uncertainty in the extrapolation. 
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Exciting new results [1| from the Sudbury Neutrino 
Observatory (SNO) have proven for the first time that 
the measured high-energy neutrino flux from the Sun 
agrees well with the one calculated from standard solar 
models 0; 01 non-electron flavour neutrinos are taken 
into account. This again focusses attention onto the 
^Be(p,7)^B reaction which provides almost exclusively 
the high-energy neutrinos measured in the SNO exper- 
iment. Their flux depends linearily on the ^Be(p,7)^B 
cross section at solar energies. Very recently, the latter 
has been redetermined by new high-precision direct mea- 
surements 0, IE IE 13 and extrapolated to zero energy 
with the help of a theoretical model 0. The resulting 
zero-energy astrophysical S factors, S'i7(0), however, do 
not always agree within their quoted errors: Hammache 
et al. 4] found 5i7(0) = 18.8 ± 1.7 eV b, in agreement 
with other direct-capture data [E IE In contrast, 

Junghans et al. 0| report a considerably larger value, 
Sn{0) = 22.3 ± 0.7 ± 0.5 eV b. The very recent resuh 
of Baby et al. also favours a rather large value of 
S'i7(0) = 21.2 ±0.7 eV b. 

In view of their importance for astro- and elementary- 
particle physics, these conflicting results should be ver- 
ified and cross-checked by other, indirect measurements 
that have different systematic errors. One possibility is 
Coulomb dissociation (CD) of ®B in the electromagnetic 
field of a high-Z nucleus. Such measurements have been 
performed at low j intermediate 0, 0| , and high en- 
ergies .14|. Alternatively, 5*17(0) can also be calculated 
from asymptotic normalization coefficients (ANC) which 
in turn are determined in low-ener g:y p roton-transfer or 
in proton-removal reactions 0, lla Il7| . 



In the present Letter, we focus on a crucial question 
that must be answered if one wants to use the CD method 
to derive a precise value for 517(0). The astrophysical S 
factors of the ''Be(p,7) reaction can only be calculated 
reliably from the energy-differential CD cross sections if 
the electromagnetic multipole components relevant for 
direct capture and the time-reversed process have the 
same strength. In low-energy proton capture the El con- 
tribution by far dominates the cross section. While El is 
the dominant multipolarity also in CD, one can show eas- 
ily that the equivalent photon field emitted from a high-Z 
target nucleus contains a strong E2 component. This is 
particularly true for CD at low energies. At higher ener- 
gies (see Ref. 0|) the relative amount of E2 multipolar- 
ity is expected to be reduced, but may still be substantial 
enough to affect the final result. To remove this ambi- 
guity, it is indispensable to either determine the E1/E2 
ratio in CD experimentally, or to extract S'17 with such 
cuts that any E2 contribution is negligible. 

Experimental limits for a possible E2 contribution were 
extracted in the work of Kikuchi et al. J^] and Iwasa et 
al. 14]. Both papers found negligible E2 contributions. 
Recently, Davids et al. have reported positive experi- 
mental evidence for a finite E2 contribution in CD of ^B, 
mainly from the analysis of inclusive longitudinal mo- 
mentum (p||) spectra of ''Be fragments measured at 44 
and 81 A MeV J^. The asymmetries in the p|[ spectra 
were interpreted to be due to E1-E2 interference in terms 
of first-order perturbation-theory . 

In order to resolve these discrepancies, we decided to 
perform an exclusive CD experiment at high energy (254 
A MeV) at the kaon spectrometer KaoS at GSI [13 with 




FIG. 1: Vector diagram showing the definitions of the an- 
gles 9cm and ^cm as well as the proton in-plane transverse 
momentum, pj", in the frame of the *B* system. 



the aim to measure quantities that should be sensitive 
to contributions of E2 multipolarity, namely the angu- 
lar correlations of the ^B-breakup particles, proton and 
^Be. Experimentally, this requires high-resolution mea- 
surements of the positions and angles of the incident ^B 
beam as well as those of the breakup fragments. The ^B 
secondary beam was produced at the SIS/FRS radioac- 
tive beam facility at GSI [23| by fragmenting a 350 A 
MeV ^^C beam in a 8 g/cm^ Be target and separating 
it from contaminant ions in a 1.4 g/cm^ wedge-shaped 
Al degrader placed in the FRS intermediate focal plane. 
Typical ^B beam intensities in front of KaoS were 5 x 10* 
per 4 sec spill; the only contaminant consisted of about 
20% ^Be ions which could be identified event by event 
with the help of a time-of- flight measurement. 

Positions and angles of the secondary beam incident 
on the Pb breakup target were measured with the help 
of two parallel-plate avalanche counters (PPAC) located 
at 308.5 cm and 71 cm upstream from the target, respec- 
tively. The detectors, which were designed at RIKEN 
plf. had areas of 10 x 10 cm^ and allowed to track the 
incident *B beam with about 90% efficiency and with 
position and angular resolutions of 1.3 mm and 1 mrad, 
respectively. Downstream from the Pb target (which con- 
sisted of 50 mg/cm^ 208pb enriched to 99.0±0.1%), the 
angles and positions as well as the energy losses of the 
outgoing particles were measured with two pairs of Si 
strip detectors (300 /im thick, 100 /im pitch) located at 
distances of about 14 cm and 31 cm. Proton and ^Be mo- 
menta were analyzed with the KaoS spectrometer which 
was set up almost identical to our previous experiment 
[T^. except for a newly constructed plastic-scintillator 
wall near the KaoS focal plane with 30 elements (each 7 
cm wide and 2 cm thick) used for trigger purposes. 

The coincident p and ^Be signals resulting from 
breakup in the ^°*Pb target were identified by recon- 
structing their vertex at the target, this removed all 
breakup events in layers of matter other than the tar- 
get. The measured momentum vectors of the outgoing p 
and ''Be particles allowed to construct the invariant-mass 
spectrum of the excited ^B* system prior to breakup. 
Fig. n shows the coordinate systems used. The angle 9^ 
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FIG. 2: In- plane transverse momenta, pi", of the break-up 
protons for three different cuts in 6s. The theoretical curves 
(full lines: El multipolarity, dashed lines: E1-I-E2 multipolar- 
ity) have been calculated in first-order perturbation theory. 
They were normalized individually to the data points in each 
frame. 



is the laboratory scattering angle of B* relative to the 
incoming ®B beam. The polar angles, 9cm, and the az- 
imuthal angles, 0cm 7 of the breakup protons are measured 
in the rest frame of the *B* system, as shown in Fig. ^ 
In the same way, one can calculate e.g. the transverse 
proton momentum vector in the reaction plane (p™). 

In the following we will present some angular distri- 
butions of the emitted proton in the frame of the ^B* 
system that can be shown to be sensitive to an E2 ampli- 
tude in CD. To interprete the measured distributions we 
need guidance by a theoretical model. To this end, we 
have performed standard first-order perturbation-theory 
(PT) calculations of the CD process in the semi-classical 
approach using a simple Woods-Saxon potential 

model for ^B. The potential depth for the ground state 
of ^B was adjusted to match the proton binding energy; 
the potential depths of the scattering states were fitted 
to the scattering lengths of the ''Li+n mirror system [23| . 
We used a radius parameter of vq = 1.25 fm and a dif- 
fuseness of a = 0.65 fm. For channel spin 1 = 2 (the 
dominant contribution) we obtained a potential depth of 
V2 = 52.60 MeV. The resulting scattering length for this 
channel of Uq'^'^" = — 8 fm agrees well with the recently 
measured value of aQ2^ = — 7 ± 3 fm (Angulo et al. j2S]). 

To take into account absorption due to nuclear overlap 
in CD, we have introduced a diffuse absorptive nuclear 
potential with a depth of 20 MeV and a radius of 9.91 
fm, i.e. the sum of the projectile and target radii. This 
choice reproduces well the integral 6*8 angular distribu- 
tion. Technically, the results of the PT calculations were 
returned as a statistically-distributed ensemble of 500 000 
CD- "events" that were analyzed in the same way as the 
experimental data, thus imposing the experimental cuts. 

We first present in Fig. [21 the distribution of for 
three different upper limits in 9s, 0.62°, 1.0°, and 2.5°. 
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FIG. 3: Top: Experimental distributions of the proton az- 
imuthal angular ((/f>cm) distributions for three different bins of 
the p-^Be relative energy, Erei- The full histograms denote a 
first-order perturbation-theory calculation for El multipolar- 
ity, the dashed ones for E1-I-E2. All theoretical curves were 
individually normalized to the data points in each frame. Bot- 
tom: the same for the polar breakup angles, 9cm- 



In classical Rutherford scattering, this corresponds to im- 
pact parameters of 30 fm, 18.5 fm, and 7 fm, respectively. 
Relative energies between p and ''Be up to 1.5 MeV were 
selected. The experimental data for all three ^g-cuts can 
be reproduced well by a PT calculation that includes only 
El multipolarity (full histograms in Fig.|21 the theoretical 
curves were normalized individually to the data points). 
If El-plus-E2 multipolarity is used in the PT calcula- 
tion, the different impact-parameter dependences of El 
and E2 multipolarity lead to markedly different shapes 
for the different 6'8-cuts (dashed histograms in Fig. 
The latter distributions are, however, in clear disagree- 
ment with our data points. 

Fig. 13 depicts the experimental 0cm and 9cm distri- 
butions for three different Erei bins, as indicated in the 
figure. A "safe" 6s limit of 1° was chosen. As expected, 
these distributions are mostly isotropic at low Erei (hi- 
dicative of s-waves) and become increasingly anisotropic 
for larger values (contributions from d- waves). For the 
(j)cm distributions, which are most sensistive to E2 admix- 
tures, the calculations for pure El multipolarity clearly 
fit best; inclusion of an E2 component shifts the maxima 
away from 90° and 270° with increasing Erei, while at 



the same time the anisotropy is reduced. Similar con- 
clusions can be drawn from the bottom part of Fig. |31 
where the proton polar angular {9cm) distributions are 
shown. The low-Erd bins show little sensitivity to E2 
multipolarity, whereas inclusion of E2 leads to a marked 
discrepancy near cos(^?cm) = 1 for the highest Erd bin. 
More detailed calculations show that at most E2 ampli- 
tudes of ^0.3 times the theoretical one from our simple 
model are simultaneously compatible with all our mea- 
sured observables. Since this would correspond to E2- 
contributions to the cross sections of less than 1%, much 
less than the errors of the data points, we neglect the ef- 
fect of E2 multipolarity. This is in line with conclusions 
drawn by Kikuchi et al. ^3 a-^d by Iwasa et al. Q from 
their respective 9^ distributions (which are, however, less 
sensitive to a small E2 component than the present angu- 
lar correlations). Our findings contradict the conclusions 
of Davids et al. W\ that a substantial E2 cross section 
has to be subtracted from the total measured CD cross 
section. 

Our results allow to interpret the relative-energy dis- 
tributions of the breakup particles in an easy way. In 
the following, we have restricted the angles 9^ to values 
below 0.62° to ensure both dominance of CD and reduc- 
tion of the effect of any possible E2 contribution. The 
data are compared to a simulation with GEANT that in- 
cludes two electromagnetic multipole components: a res- 
onant Ml contribution located at £'rei=0.63 MeV with 
resonance parameters taken from Filippone et al. |lOj |. 
and the non-resonant El contribution from our theoreti- 
cal model as described above. The latter was scaled by a 
normalization factor of 0.79. Note that we have added to 
the GEANT simulation a contribution that feeds the first 
excited state in "^Be at 429 keV using the measurements 
of Kikuchi et al. Subtracting the small Ml contri- 

bution (that affects only a narrow Erei region around the 
resonance), the remaining da/dErei distribution can be 
converted to the El astrophysical S factor Sn{Erei). 

The resulting 6*17 factors (averaged over Erei bins 0.2 
to 0.3 MeV wide) are visualized in Fig. 0] The error 
bars do not include a common systematic error of 5.6%. 
The top panel (a) compares our results to those of other 
CD experiments [H 111 111 (the data of Ref.jil rep- 
resent their F,l-Sn factors after subtraction of the E2- 
contribution) . At low Erei , the CD S factors are in good 
agreement, though the Davids et al. ^3 data are sys- 
tematically lower. The bottom panel (b) compares our 
data to those of the ''Be(p,7)®B measurements where the 
authors have subtracted the contribution from the Ml 
resonance (Refs. j^i^^]). At low energies the (p,7) data 
of Refs. 0, 0] and ours are in good agreement, whereas 
the Seattle data @ deviate considerably. The opposite 
behaviour is noted above the Ml resonance: our data 
and those of Refs. 0, Q match excellently, whereas the 
other (p,7) experiments 0, 0,0, 0| consistently report 
lower values. We want to emphasize the remarkably good 
agreement of our CD data up to 1.1 MeV with the most 
recent direct-proton-capture experiment where an ion- 



4 



> 



30 



20 



10 
40 



30 



20 



10 



• this work (GSI-2) 

O Iwasa et al. 

M Kikuchi et al. 

□ Davids et al. (E1 only) 




b) 



n 1 r 

• this work (GSI-2) 



□ Hammache et al. 
* Junghans et al. 
O Baby et al. 



this work 




Descouvemont 



0.25 



0.5 



0.75 



1.25 



1.5 



FIG. 4: a) Comparison between Sn values from Coulomb- 
dissociation experiments. The full (open) circles indicate the 
present (previous) GSI CD experiment. Open stars depict 
Ref. [l^. open squares Ref. (E2 contribution subtracted). 
The theoretical curves are described in the text, 
b) 5i7 from this work in comparison with the (p, 7) experi- 
ments of Ref . Q (squares), Ref. (stars), and Ref. j3 (open 
circles). The latter data were corrected for the contribution 
of the Ml resonance by the authors. 



implanted '^Be target was used 0- 

To extrapolate to zero energy, all recent (p,7) exper- 
iments have chosen the cluster model of Descouvemont 
and Baye . When we fit our data points up to Erei — 
1.5 MeV to this model and add in quadrature a common 
systematic error of 5.6%, we obtain S'i7(0) = 20.8 ±1.3 
eV b (dashed lines in Fig. 0}. Restricting the fit to en- 
ergies below 0.6 MeV, where the model-dependence has 
been shown to be weaker [H, £'17(0) = 19.6 ± 1.4 eV 
b is obtained. Our potential model, however, repro- 
duces the data over the entire energy range up to 1.5 
MeV, yielding 517(0) = 18.6 ± 1.2 eV b (full hues in 
Fig. 0J . It is interesting to note that a fit of the Baby et 
at (p,7) data to our model yields practically the same 
result, 517(0) = 18.1 ± 0.3 eV b. Clearly, still more 
high-precision experimental data are needed to resolve 
the discrepancies between the experimental data sets and 
to pin down the correct theoretical extrapolation of the 
measured data to solar energy. In the meantime, an ad- 
ditional "extrapolation error" of ± 1.0 eV b seems ap- 
propriate. 

We conclude that Coulomb dissociation has been 
proven to be a valuable method to provide a rather pre- 
cise value for the low-energy '^Be(p,7) cross section. Since 
in CD all energy bins are measured simultaneously, CD 
provides a reliable measurement of the shape of the Sn 
distribution. By setting tight constraints to the scatter- 
ing angle 6^ and analyzing p-^Be angular correlations, 
a significant contribution from E2 multipolarity can be 
excluded. Small modifications of the Woods-Saxon po- 
tential parameters allow to reproduce the data in first- 
order perturbation theory with remarkable accuracy up 
to about Erei = 1.5 MeV. 
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